Equilibrium shapes of flat knots 
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We study the equilibrium shapes of prime and composite knots confined to two dimensions. Using 
rigorous scaling arguments we show that, due to self-avoiding effects, the topological details of prime 
knots are localised on a small portion of the larger ring polymer. Within this region, the original 
knot configuration can assume a hierarchy of contracted shapes, the dominating one given by just 
one small loop. This hierarchy is investigated in detail for the fiat trefoil knot, and corroborated by 
Monte Carlo simulations. 
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The static and dynamic behaviour of single polymer 
chains, such as DNA, and multichain systems like gels 
and rubbers, is strongly influenced by knots and perma- 
nent entanglements |[ ||, ^. Topological constraints 
are created with probability one during the polymerisa- 
tion of long closed chains (5); more generally, knots and 
entanglements are a ubiquitous element of higher molec- 
ular multi-chain melts and solutions. This has profound 
consequences, reaching far into biology and chemistry. 
For instance, knots in DNA impede the separation of 
the two strands of the double helix during transcription, 
and therefore the access to the genetic code Q . Chemi- 
cally, even single closed polymers may exhibit quite dif- 
ferent properties if they have different topology ||^ . In the 
nanosciences, recent experimental techniques allow sin- 
gle polymer molecules (with fixed topology) to be probed 
and manipulated; e.g., by single molecule spectroscopy 
1^ , or by optomicroscopical imaging of small latex beads 
attached to a molecule These tools provide impe- 
tus for the theoretical understanding of the behaviour 
of macromolecules under topological constraints. How- 
ever, analytical studies, such as the statistical mechan- 
ics of a knotted polymer, are difficult since topological 
constraints require knowledge of the complete shape of 
the curve. Such global constraints are hard to imple- 
ment, and a complete statistical mechanical description 
of knots remains unattained ^ ^ . 

The mathematical discipline of knot theory provides 
tools for the classification of knots. In particular, differ- 
ent knots can be distinguished by their projections onto 
a 2D plane, keeping track of crossings according to which 
segment passes on top of another Q . By a sequence of 
so-called Rcidcmcister moves [||, which leave the topol- 
ogy unchanged, the number of crossings can be reduced 
to a minimum, which is a simple topological invariant 
1^, For instance, in Fig. ^ we depict the minimal pro- 
jection of the trefoil knot, classified as 3i, with its 3 cross- 



ings. Such quasi 2D projections, which we call flat knots, 
can be physically realized by compressing originally 3D 
knots by forces normal to the projection plane. Examples 
include polymers adsorbed on a surface or membrane by 
electrostatic or other adhesive forces [|3| ; or confined be- 
tween parallel walls. In these cases the flat polymer knot 
can still equilibrate in 2D. Another experimental realiza- 
tion comes from Ref. ]l^ , in which macroscopic knotted 
chains are flattened by gravity onto a vibrating plane. 
The equilibrium shapes, and their scaling properties, of 
such flat knots are studied in this paper. Flat knots have 
the additional advantage of being easy to image by mi- 
croscopy. They are also more amenable to numeric stud- 
ies than their 3D counterparts, and have in fact been 
already studied in Ref. jl^ . 

There is growing numerical evidence that prime knots 
are tight in the sense that the topologically entangled re- 
gion is statistically likely to be localised on a small por- 
tion of the longer chain jl^, |l6|, . Indirect numerical 
evidence of this was originally obtained by simulations in- 
dicating that the radius of gyration of a long polygon in 
3D is asymptotically independent of its knot type; while 
the presence of the knot increases the number of config- 
urations by a factor related to the number of positions of 
the tight region around the remaining loop fl^ . Simula- 
tions of 2D polygons in Ref. jist provide quite convincing 
visual evidence of localised knot regions. In this paper, 
we quantify the tightness of flat knots, using scaling ar- 
guments to obtain the power law size distributions for a 
hierarchy of possible equilibrium shapes. For the trefoil. 
Fig. shows this hierarchy of shapes and the correspond- 
ing exponents for the distribution of knot size. 

To demonstrate the entropic origin of tight shapes, we 
initially consider a simple ring of length L with one cross- 
ing. The effect of the crossing consists in creating two 
loops of lengths £ and L — £, respectively, while the ori- 
entation of the crossing is irrelevant. In this sense, the 
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Figure 1: Standard minimal projection of the trefoil knot 3i, followed by its different possible contractions, arranged according 
to higher scaling orders. The uncontracted trefoil geometry is found at position III of the hierarchy. At I, the figure-eight 
structure is drawn. The diagrams II- VI show the multiply connected knot region of total length £ = Ylii^^ where the 
protruding legs indicate the outgoing large loop of length sjk/ = L — £. Below the individual contractions, we include the 
network exponents 75, the number m of independent integrations, and the exponents c defined via the PDF, p{£) £~'^. 



crossing can be replaced with a vertex with four outgoing 
legs. We denote the resulting network by Qi (see posi- 
tion I in Fig.|l|). In fact, we can more generally consider 
a sliding ring, or slip- link to force two points of the 
chain to be close to each other, in d dimensions. With- 
out self-avoiding constraints (ideal chains), the number 
of configurations lui{£, L) scales as |l9| 

LOi{£,L)^ ^l^£-'^'^{L-£)-'^'^ , (1) 

where, on a lattice, /i is the effective connectivity con- 
stant for Gaussian random walks. The average loop 
size is given by {£) = J^'" d££uji{i, L) j /^^"° dluji{£, L), 

where a is a short-distance cutoff set by the lattice con- 
stant. Note that {£) = L/2 due to symmetry. However, 
the corresponding probability density function (PDF) is 
strongly peaked at ^ = and i = L, and a typical shape 
consists of one tight and one large loop. In = 2, the 
mean size of the smaller loop, (£)< ~ L/\\n{a/L)\, is 
still rather large. It is instructive to compare to higher 
dimensions: one has weak localisation, (£)< ^ a^/^L^^^, 
in d = 3, and strong localisation, (^)< ^ a, m d > A. 
Thus, for ideal chains, tightness of the smaller loop is 
more pronounced in higher dimensions. 

To include self-avoiding interactions, we use results for 
general polymer networks obtained by Duplantier [ pO[ , 
and in Refs. [|l| In a network Q consisting of M 
chain segments of lengths si , . . . , sj\f and total length L = 
X]i=i number of configurations Lug scales as 

u;g{s,, ...,sm)^ l^'-sl^-'yG (—,■■■, , (2) 

where yg is a scaling function, and /i is the effective con- 
nectivity constant for self-avoiding walks. The exponent 
^g is given by 75 = 1 — dvC + X! ^n<Jn, where v is 

N>1 

the swelling exponent, C is the number of independent 
loops, riAT is the number of vertices with N outgoing legs. 



and (Tat is an exponent associated with such a vertex. In 
d = 2, (TAT = (2 - N){m + 2)/64 0. 

The network Qi corresponds to the parameters Af — 2, 
C — 2, — I, si — £, and S2 = L — £. By virtue of 
Eq. (|^), the number of configurations of Qi with fixed t 
follows the scaling form 

u;,{£,L)^^,\L~£r^-'x(j^^, (3) 

where 71 = 1 — 2dv + 0-4. In the limit i <^ L, uJi{£,L) 
should reduce to the number Wcrw(i) ~ iJ,^L~'^'^ of closed 
random walks of length L which start and end at a given 
point in space ||l9|, This implies X{x) ~ jjTi-i+d'' ^s 
X 0, such that 

u;i{£,L)^ fi^{L~£)-'^''£-\ £ <^ L , (4) 

where c = —(71 — 1 + dv) = dv — (74. Using (74 ~ —19/16 
and u = 3/4 in d = 2, we find c = 43/16 = 2.6875. 
In d = 3, (74 « -0.48 and v « 0.588, so that c « 2.24 
1^ , |2^ . In both cases the result c > 2 imphes that the 
loop of length £ is strongly localised in the sense defined 
above. This justifies the a priori assumption £ L, and 
makes the analysis self-consistent. Note that for self- 
avoiding chains, in d = 2 the localisation is even stronger 
than in d = 3, in contrast to the corresponding trend for 
ideal chains. 

We performed Monte Carlo (MC) simulations (see be- 
low for details) of a chain in d = 2 confined to the Qi 
(figure-eight) structure by a slip-link. As shown in Fig. 
^, the size distribution for the small loop can be fitted 
to a power law with exponent c = 2.7 ± 0.1 in good 
agreement to the above prediction. 

For the configuration Qj, the probability for the size 
of each loop is peaked a,t £ ^ and £ — > L. For more 
complicated projections, the joint probability to find the 
individual segments with given lengths Si is expected to 
peak at the edges of the higher-dimensional configuration 
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Figure 2: Power law tails in PDFs for the size i of tight 
segments: As defined in the figure, we show results for the 
smaller loop in a figure-eight structure, the overall size of 
the trefoil knot, as well as the two leading contractions of the 
latter. The insets show typical configurations of the small loop 
for a figure-eight (the arrow points to the slip-link consisting 
of 3 tethered beads) , and the trefoil. 

hyperspace. Some analysis is necessary to find the opti- 
mal shapes; as presented here for the simplest non-trivial 
knot, the (flat) trefoil knot 3i (see Fig.|l|). Each of the 
three crossings is replaced with a vertex with four outgo- 
ing legs, and the resulting network is assumed to separate 
into a large loop and a multiply connected region which 
includes the vertices. Let I = Y^i=\ be the total length 
of all segments contained in the multiply connected knot 
region (see Fig.|[ position III). Accordingly, the length 
of the large loop is L — £. In the limit £ -C L, the number 
of configurations of the network ^ni can be derived in a 
similar way as above, yielding 

c.J„ ^ M^(L - ^)-^^r--i+'^''>V (^, |, |, ^) , (5) 

where 7111 = 1 — /^dv + 3(74 and W is a scaling function. 
The prime on wm indicates that each of the segment 
lengths Si is fixed. In order to obtain the number of 
configurations a;iii(£, L) for the case we are interested in, 
where only the total length t is fixed, we should inte- 
grate lo'yyi over all distributions of lengths Si under the 
constraint X)i=i — ^- This leads to the result 

u;III(£,L)~A^^(L-^)-'*''r^ (6) 

with c = — (7111 — 1 -l- dv) — m, where m = 4 corresponds 
to the number of independent integrations over Si. Thus, 
c = Sdi' — 3f74 — 4 = yI (see Fig.|^, position III). 

However, some care is necessary in performing these 
integrations, since the scaling function W in Eq. (^ may 
exhibit non-integrable singularities if one or more of its 
arguments tend to or 1. The geometries corresponding 



to these limits (edges of the configuration hyperspace) 
represent contractions of the original trefoil network Qm 
in the sense that the length of one or more of the seg- 
ments Si is of the order of the short-distance cutoff a. 
If such a short segment connects different vertices, they 
cannot be resolved on larger length scales, but melt into 
a single, new vertex, in the context of our scaling analy- 
sis . Thus, each contraction corresponds to a different 
network Q, which may contain a vertex with six or even 
eight outgoing legs. For each of these networks, one can 
calculate the corresponding exponent c in a similar way as 
above, and using the Euler relations 2A/' = J2n>i 
and C = X]jv>i ^ 2)njv + 1, we obtain 

ce = 2 + En>4 [§{diy - 1) + {\aN\ - di^)] . (7) 

Our scaling analysis relies on an expansion m a/£ <^ 1, 
and the values of c determine a sequence of contractions 
according to higher orders in a/£: The smallest value of 
c corresponds to the most likely contraction, while the 
others represent corrections to this leading scaling be- 
haviour, and are thus less and less probable (see Fig. |^). 
To lowest order, the trefoil behaves like a large ring poly- 
mer at whose fringe the point-like knot region is located. 
At the next level of resolution, it appears contracted to 
the figure-eight shape (see Fig. 0, position I) . If more ac- 
curate data are available, the higher order shapes II to 
VII may be found with decreasing probability. Inter- 
estingly, the original uncontracted trefoil configuration 
ranks third in the hierarchy of shapes. Note that the 
contractions shown in Fig. |l] may occur in different topo- 
logical variants. For instance, the smaller loop in contrac- 
tion I could be inside the larger loop. However, this does 
not make a difference in terms of the scaling analysis. 

To check these predictions, we performed MC simu- 
lations of a flat trefoil knot using a standard bead-and- 
tether chain model in 3D. The polymer was flattened by 
a "gravitational" fleld V = —ksTh/h* perpendicular to 
a hard wall, where h is the height and h* was set to 0.3 
times the bead diameter. With 512 beads, around 10^^ 
MC steps were performed to generate equilibrium states. 
Configurations corresponding to contraction I are then 
selected by requiring that besides a large loop, they con- 
tain only one segment larger than a preset cutoff length 
(taken to be 5 monomers), and similarly for contraction 
II. The size distributions for such contractions, as well 
as for all possible knot shapes are shown in Fig. ^. The 
tails of the distributions are indeed consistent with the 
predicted power laws, although the data (especially for 
contraction II) is too noisy for a definitive statement. 

These scaling results pertain also to other prime knots. 
In particular, the dominating contribution for any prime 
knot corresponds to the figure-eight contraction Qi. This 
can be shown by noting that Eq. (|^) predicts a larger 
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Figure 3: Typical equilibrium configurations for the trefoil 3i, the prime knot 819, and the composite knot 3i#3i consisting of 
two trefoils, in d = 2. The initial conditions were symmetric. 



value of the scaling exponent c for any network Q other 
than Qi . Figure ^ demonstrates the tightness of the prime 
knot 819. Composite knots, however, can maximise the 
number of configurations by splitting into their prime fac- 
tors as indicated in Fig.^ for 3i#3i. Each prime factor 
is tight and located at the fringe of one large loop, and 
accounts for an additional factor of L for the number of 
configurations, as compared to a ring of length L without 
a knot. Indeed, this gain in entropy leads to the tightness 
of knots. 

In conclusion, we find that the trefoil knot, as well as 
higher order prime and composite knots, are sharply lo- 
calized when forced to lie flat. In the most likely shapes, 
each prime factor is tightened into a loop (a figure-eight 
contraction) . It is natural to speculate that entropic fac- 
tors also confine knots in d = 3 by power law distributions 
in size. Direct checks of such behaviour are hampered 
by the difficulty of identifying the knotted region of a 
curve in d = 3. One may instead search for indirect 
signatures of localized knots in detailed dependencies of 
gyration radius and other polymeric quantities on length. 
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